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Abstract. Exclusive p'^-meson electroproduction is studied in the HERMES experiment, using a 27.6 GeV 
longitudinally polarized electron/positron beam and unpolarized hydrogen and deuterium targets in the 
kinematic region 0.5 GeV^ < < 7.0 GeV^ 3.0 GeV < < 6.3 GeV, and -t' < 0.4 GeV^ Real 
and imaginary parts of the ratios of the natural-parity-exchange helicity amplitudes Tn (7^^ — >• p^), Toi 
ilT ~^ Pl), Tio {-yl — >■ Pt), and Ti-i {'jIt ~^ Pt) to Too (7l ~^ Pl) are extracted from the data. For the 
unnatural-parity-exchange amplitude Un, the ratio |?7ii/Too| is obtained. The and t' dependences of 
these ratios are presented and compared with perturbative QCD predictions. 
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1 Introduction 

Exclusive electroproduction of vector mesons, e + — >■ 
e' +V + N' , has been the focus of investigation for decades. 
Not only is the reaction mechanism of intrinsic interest, 
but this process also offers the possibility of studying, in a 
model-dependent way, the structure of hadrons involved in 
the process [TJ[2] . Using the one-photon-exchange approx- 
imation, all the measurable observables in electroproduc- 
tion can be expressed in terms of the virtual photon spin- 
density matrix and the helicity amplitudes ^AvA'„a.^Ajv of 
the process 

r{X^) + N{XN)^V{Xv) + N'{yj^) . (1) 

Here, 7* denotes the virtual photon exchanged between 
the lepton and the target nucleon, V denotes the pro- 
duced vector meson, and N{N') the initial (final) nucleon. 
The hclicitics of the particles are given in parentheses in 
Eq. (Ill) . The helicity amplitudes are defined in the virtual- 
photon-nucleon center-of-mass (CM) system. In order to 
make the discussion of the transition ^ V more trans- 
parent, we shall henceforth often omit the nucleon helicity 
indices using the notation F\^,a^ instead of i^Ai^A'^A^Aw • 
The properties of helicity amplitudes can be studied in 
detail because the spin-density matrix of the virtual pho- 
ton is well known from quantum electrodynamics and the 
spin-density matrix of the produced vector meson is ex- 
perimentally accessible. 

For unpolarized targets, the formalism of the spin- 
density matrix elements (SDMEs) of the produced vector 
meson was first presented in Ref . [3] , where expressions of 
SDMEs in terms of helicity amplitudes were established. 
The formalism was then extended to the case of polarized 
targets in Ref. [J. Recently, a new general formalism for 
the description of the process in Eq. ^ through SDMEs 
was presented in Ref. |S]. 

In order to determine SDMEs from experimental data, 
the SDMEs are considered as being independent free pa- 
rameters in the fitting of the production and decay angular 
distribution of the vector meson. This is referred to as the 
"SDME method" in the rest of this work. 

SDMEs arc dimcnsionlcss quantities and therefore de- 
pend on ratios of amplitudes rather than on amplitudes 
themselves. The exact expressions for SDMEs given in 
Refs. I31IHI5] can be rewritten in terms of amplitude ratios. 



In an alternative method of fitting the angular distribu- 
tion, these ratios are considered as being independent free 
parameters. This method is referred to as the "amplitude 
method" in the rest of this work. 

In order to extract the helicity amplitudes themselves, 
experimental data on the differential cross section with re- 
spect to the Mandelstam t variable, da/dt (which is pro- 
portional to the sum of squared moduli of all the ampli- 
tudes) are required in addition to the experimentally de- 
termined amplitude ratios. An analysis of these combined 
data would allow the extraction of the moduli of all am- 
plitudes and of the phase differences between them with 
the common phase remaining undetermined. However, the 
requisite information on nucleon spin-flip amplitudes is 
not available for data taken with unpolarized targets. 

Exclusive meson production in hard lepton scattering 
has been shown to offer the possibility of constraining gen- 
eralized parton distributions (GPDs), which provide cor- 
related information on transverse spatial and longitudinal 
momentum distributions of partons in the nucleon |S1[71 
liinilTnEIlIIllIllIlIISlIIS- Vector-meson production am- 
plitudes contain various linear combinations of process- 
independent GPDs for quarks of various flavors and glu- 
ons. Access to GPDs relies on the factorization property 
of the process amplitude, i.e., the amplitude can be writ- 
ten as convolution of "non-perturbative" GPDs with am- 
plitudes of hard partonic subprocesses calculated within 
the framework of perturbative quantum chromodynamics 
(pQCD) and quantum electrodynamics. 

The amplitudes -FqIq-i-i are the most interesting be- 
cause, only for these and quantities constructed from them, 
was factorization proven |17j . Factorization is not proven 
for the other amplitudes nor for SDMEs, which depend 
on all helicity amplitudes. The amplitudes -Fgig^i corre- 
spond to the transition of a longitudinally polarized (L) 
virtual photon to a longitudinally polarized vector meson, 
7l ~^ ^L, and dominate at large photon virtuality Q^. 
Nevertheless, an application of the "modified perturbative 
approach" [T^lfn] assumes that the factorization property 
also holds for the amplitudes Fu and Fqi. The agreement 
found between certain calculated SDMEs and those ex- 
tracted from HERMES [IS], ZEUS 01, and HI gD] data 
supports this assumption. The differential and total cross 
sections for longitudinal p^-meson production by longitu- 
dinal photons are reasonably well described in the GPD- 
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based approach of Refs. [ni[T^[TIl[T51[TBl not only at the 
high energies of the HERA colhder experiments [2n[22l 
[23, 24,25] but also at intermediate energies covered by the 
fixed-target experiments E665 g^, HERMES [37], and 

CLAS [MUH]. 

Amplitude ratios are useful not only as a tool for ob- 
taining Fqq. They can be used to study the general prop- 
erties [51 I30[I5T1I5^ of the vector-meson-production ampli- 
tudes at very small and very large values of and t. 
They can also be used to test theoretical models. For these 
purposes, an amplitude ratio is more convenient than an 
SDME, which constitutes a more complicated object since 
any SDME depends on all the amplitude ratios. The am- 
plitudes of the non-diagonal V transitions Fqi and 
Fio are known |31j to be zero if any valence quark in 
the vector meson carries half of the longitudinal momen- 
tum of the meson in the infinite momentum frame. These 
amplitudes are the most sensitive objects for the study 
of the quark motion in vector mesons. The double-spin- 
flip amplitude i^i_i contains information on the gluon 
transversity generalized parton distribution [Ql lTOllSOllSTl 
[531155] in the nucleon, which cannot be obtained from in- 
clusive deep-inelastic Icpton-nucleon scattering. 

Amplitude ratios in p° and (j) meson production on 
the proton were first studied by the HI experiment [20] 
at the HERA collider. The results of the analysis of 
meson production at HERMES using the SDME method 
were published in Ref. [TH], where also the description of 
the HERMES spectrometer and details of the data treat- 
ment etc. can be found. The present work is a continuation 
of that analysis. In the analysis presented in this paper, 
ratios of hclicity amplitudes to the amplitude ^0^0^ 
extracted separately for data taken with unpolarized hy- 
drogen and deuterium targets first and then also for the 
combined datasets in a single analysis. A comparison of 
the proton and deuteron results allows the study of the 
degree of interference between 1=1 exchanges of qq pairs 
and 1 = exchanges of qq pairs and two gluons. 



2 Kinematics 

In accordance with the notation of Ref. [TH] , the kinematic 
variables of the process 

eN e'7*iV ^ e'p°N' e'7r+TT~N' 



(2) 



are defined as follows. The four-momenta of the incident 
and outgoing leptons are denoted by k and A:', the differ- 
ence of which defines the four-momentum g = fc — fc' of 
the virtual photon 7*. In the target rest frame (which is 
also called the laboratory or lab frame in this work), 1) 
is the scattering angle between the incident and outgoing 
leptons, the energies of the leptons are denoted by E and 
E' . The photon virtuality, given by 



-q^ = -{k-k 



,^ o lab , o 17 

'^2 ^ AEE'sin^ 



2' 



(3) 



is positive in leptoproduction. In this equation, the elec- 
tron rest mass is neglected. The four-momentum of the in- 
cident (recoiling) nucleon is denoted by p (p'). The Bjorken 
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Fig. 1. Definition of angles in the process eA*" — > e'^*N — >■ 
e' p°N' — >■ eV^Tr" A'^'. Here $ is the angle between the pro- 
duction plane and the lepton scattering plane in the CM system 
of virtual photon and target nucleon. See Ref. |34) for details. 
The variables 6 and (f) are respectively the polar and azimuthal 
angles of the decay tt"'" in the p°-meson rest frame, with the Z 
axis being anti-parallel to the outgoing nucleon momentum. 



scaling variable xb is defined as 



xb 



with 



2p-q 2Mv' 



P ■ 1 lab p 

M 



(4) 



(5) 



so that v represents the energy transfer from the incoming 
lepton to the virtual photon in the laboratory frame. The 
squared invariant mass of the photon-nucleon system is 
given by 



{p + qf = + 2Mv-Q'^. (6) 
The Mandelstam variable t is defined by the relation 

t={q-v)\ (7) 

where v is the four- momentum of the p*^ meson being equal 
to + Pt^- , the sum of 7r+ and tt^ four-momenta. The 
variables to, and t' = t — to are always negative, where 
—to is the minimal value of — t for given values of Q^, W, 
and the p°-meson mass My- At small values of ~t' , the 
approximation —t'^iv^ holds, where vt is the transverse 
momentum of the meson with respect to the direction 
of the virtual photon in the CM system. 

The variable e represents the ratio of fluxes of longitu- 
dinally and transversely polarized virtual photons and is 
given by 



1-2/ 



lab 



1 + 2 1 



tan 



(8) 
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with y ~ p ■ q/p ■ k '= h'/E. 

The "exclusivity" of p° production in the process in 
Eq. (0) is characterized by the variable 



AE- 



Ml 



lab 



2M 



Ey - [E^^ 



E„ 



(9) 



where Mx = \/{k — k' + p — p^+ — p^- )^ is the recon- 
structed invariant mass of the undetected hadronic system 
(missing mass), Ey = v + t/{2M) is the energy of the ex- 
clusively produced p° meson, and {E^+ +E^~) is the sum 
of the energies of the two detected pions in the laboratory 
system. For exclusive p*^ clectroproduction Mx ~ M and 
therefore AE = 0. 

The angles used for the description of the process in 
Eq. (0) are defined according to Ref. [M] (see also Ref. [T5]) 
and presented in Fig. [TJ 



3 Formalism 

3.1 Natural and Unnatural-Parity Exchange Helicity 
Amplitudes 

In Ref. [3] , exclusive p°-meson Icptoproduction is described 
by helicity amplitudes ^Av-a'^a^An defined in the right- 
handed CM system of the virtual photon and target nu- 
clcon. In this system, the Z-axis is aligned along the vir- 
tual photon three-momentum q and the F-axis is parallel 
to q X V where v is the p° meson three-momentum as 
shown in Fig. [TJ The helicity amplitude can be expressed 
as the scalar product of the matrix element of the electro- 
magnetic current vector J** and the virtual-photon polar- 
ization vector el^"'"' 

i^A.AWA-,A„ = {~l)^-'{vXvp'X'N\J''\P^N)ei^^\ (10) 

where a summation over the Lorentz index k is performed. 
Here el^^'* and ejf' indicate respectively transverse and 
longitudinal polarization of the virtual photon in the CM 
system: 

e(±i) = ieo,ex, eY,ez) = (0,T^,-^,0) , (11) 



=(0) _ 



^(^z, 0,0, go) 



= -^(A/v/^^^TQ2,0,0,A/t^-g2) . (12) 

The ket vector \pXn) corresponds to the initial nucleon 
and the bra vector {vXyp' X'n\ represents the final state 
consisting of a meson and the scattered nucleon. 

Any helicity amplitude -FavA'„a.^Ajv decomposed 
into the sum of an amplitude Tx^x' a^An ^^^r natural- 
parity exchange (NPE) and an amplitude Uxv\'^\-,\n ^^r 
unnatural-parity exchange (UPE) [31lll[5] 



which obey the following symmetry relations 

AvAjyA-yAjv V / — AvAjy— A^AjY 

= (-l)'~^^"T^A.-A' A.-A„, (14) 



C^Av-A'„A^An --(-1) ^''^^^^^-Ai.A'„-A.,A, 
= -(-l)''"~^"f/A.-A'„A,-A„ 



(15) 



There are three important consequences of the symmetry 
relations (HH) and ^ 

i) The number of linearly independent NPE amplitudes 
is equal to 10 while only 8 independent UPE amplitudes 
describe the process in Eq. ([T]), 

ii) No UPE amplitude exists for the transition 7^ 
so that in particular F^iqi = 



00, 



iii) For unpolarized targets there is no interference be- 
tween NPE and UPE amplitudes [3','5]. 

In Regge phenomenology [,35 .,36] . an exchange of a 
single natural- parity reggeon [P = (— 1)"^ (pomeron, sec- 
ondary reggeons p, /2, 02, ...)] contributes to the NPE am- 
plitudes, while an exchange of a single unnatural-parity 
reggeon [P = —{—lY (tt, oi, 61,...)] contributes to the 
UPE amplitudes [37] • It is worth noting that for a multi- 
reggeon-exchange contribution, there is no such one-to-one 
correspondence. For example, an exchange of two reggeons 
of "unnatural" parity contributes to the NPE amplitudes. 

For convenience, we introduced in Ref. [18) the abbre- 
viation ^ = ^X^A' Ajv ^"^^ ^^^^ summation over the final 
nucleon helicity indices and for averaging over the initial 
spin states of the nucleon. For NPE amplitudes, transi- 
tions that are diagonal in nucleon helicity (A^ = Aat) 
are dominant. In this case, neglecting the small nucleon 
helicity-flip amplitudes 1 and using Eq. ([T4[) . 

the summation and averaging operation ^ reduces to one 
term: 



^Ta,,a,T;,,a; = \^t^- 



Xy A-j 



iT* 



Txv\~,Tx'x' 



(16) 



where y y represents the complex conjugate of the 
amplitude Ta'^^a^a^An ■ to symmetry properties (see 
Eq. ([13])), the abbreviated notation TavA^ is used for the 
amplitudes T^, 1 x 1 = Tx,, 1 1, 1 . 

In general, no dominance for UPE amplitudes can be 
proven either for diagonal transitions {Xn = X'x) or for 
transitions with nucleon helicity flip. Therefore, no rela- 
tion analogous to Eq. (fTSI) can be derived. 



3.2 Spin Density Matrix Elements 

The photon-spin-density matrix normalized to unit flux of 
transversely polarized virtual photons embodies the ma- 



FxyX'^x^x^ = Ta^,a^a.,a« + Ux^x'^x^x,^, (13) trices g^^^^ and g^_^^_^ corresponding to unpolarized (U) 
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and longitudinally (L) polarized lepton beams: 



J7+L 



Pb q\ 



(17) 



where Pb is the longitudinal polarization of the lepton 
beam and and are defined in Ref. [18]. The 

spin density matrix p\vnv of the produced meson is 



related to that of the virtual photon, g\^^ , through the 
von Neumann formula [3] : 



JJ+L p* /-IQ-l 

tfX^fj.^ ^/^v A^/i^ Ajv ■ \ 



Note that the phase factor (— 1)'^^' in Eq. pUj) is important 
for usual matrix summation over and py in Eq. ()18p . 
The normalization factor is given by 



M =Nt + eJ^L 



(19) 



with 



E 



|C/ii|' + |C/oiP + |C/i_i|2 



2|C/i, 



(20) 



(21) 



Equation (|2ip is obtained by using the symmetry rela- 
tions (mi) and (dSl). 

If the spin-density matrix of the virtual photon is de- 
composed into the standard set (sec Ref. [3]) of nine her- 
mitian matrices {rj = 0, 1, 8), for the produced p*^ 
meson a set of nine matrices , is obtained: 

— ^ 'S^ F I F* ('29'! 

XvP-v 2J\f Z^^vX^XyXN ^X.jfi-,-'^ tivX'j^fi-,XN ' V / 

A^/^-y Ajy A/v^ 

If an experiment is performed at fixed beam energy, a 
(Roscnbluth |38| ) decomposition into contributions from 
longitudinally and transversely polarized virtual photons 
is impossible. In this case, the contributions from the ma- 
trix elements r? , and rf , cannot be disentangled 
and the only measurable combination is 



Xy fJ'V — fJ'V 



AvMv 



(23) 



4 The Amplitude Method 

4.1 Comparison of the SDME and Amplitude IVIethods 

The SDME method for analyzing experimental data has 
the disadvantage that, in general, numerators of SDMEs 
depend on several ratios of helicity amplitudes (see Refs. [3l 
[T5] ) while the denominator common for all SDMEs de- 
pends on all the amplitude ratios. Therefore, if an SMDE 
extracted from data differs from model calculations, the 
source of the discrepancy is difficult to identify. 

In total there are 18 independent amplitudes OHUS]. 
Since SDMEs depend on ratios of these complex ampli- 
tudes to Fqo, any SDME comprises 34 real, independent 
functions. However, the number of SDMEs to be extracted 
from experimental data in polarized particle scattering is 
larger. For instance, there are 47 SDMEs when both beam 
and target are longitudinally polarized, while the angular 
distribution for longitudinal beam polarization and trans- 
verse target polarization depends on 71 SDMEs [S|. While 
the SDMEs can be completely expressed in terms of am- 
plitude ratios, the fact that there are many more SDMEs 
than amplitude ratios implies that the SDMEs obey some 
relations. Due to the complicated inter-dependence on am- 
plitude ratios, SDMEs cannot be considered as indepen- 
dent quantities and amplitude ratios can provide a more 
economical basis for fitting angular distributions of decay 
particles. This implies that the SDME values calculated 
from extracted amplitude ratios can be expected to be 
more accurate than those from direct SDME fits. 

SDMEs calculated using the extracted helicity ampli- 
tude ratios may differ from those obtained with the SDME 
method. The SDMEs calculated from amplitude ratios 
are more constrained than those obtained by the SDME 
method. In general, any set of SDMEs obtained with the 
SDME method is physical only if it can be described in 
terms of amplitude ratios. Certain conditions, called pos- 
itivity constraints, that must be satisfied by SDMEs in 
order to be expressed through amplitude ratios were con- 
sidered in Ref. [5]. However, the full set of conditions that 
must be satisfied by SDMEs in order to be expressible in 
terms of helicity amplitude ratios is currently unknown. 



3.3 Three-dimensional Angular Distribution 

The formula for the angular distribution of the scattered 
electrons/positrons and tt^ mesons from the decay p'^ — >■ 
TT+TT" for the case of a longitudinally polarized beam and 
an unpolarized target [5] is given by 



W^+^i-P, (p, cos 9) 
W^(<?, (I), cos 6) + PBW^{^,(j}, cose 



(24) 



Here {<l>,(f>,cosd) represents the angular distribution 
when both beam and target are unpolarized (see Eq. ((25)) ) 
and yV-^(<?, (j), cos 9) is the additional term arising for lon- 
gitudinally polarized beam (see Eq. (j26p '). 



4.2 Hierarchy of Helicity Amplitudes 

The number of amplitude ratios to be extracted from data 
can be reduced if there exists a hierarchy for the moduli 
of helicity amplitudes. For large photon virtuality and 
small such a hierarchy among helicity amplitudes was 
predicted theoretically |30[I31] and confirmed experimen- 
tally [MllIHllSn] ■ According to this hierarchy, there are a 
few significant amplitudes while the contributions of other 
amplitudes may be neglected within the present experi- 
mental accuracy. In a given kinematic region, the hierar- 
chy of the amplitudes is governed by one or more "small 
kinematic parameters" that determine the contribution of 
various amplitudes to the process. 
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W^(<P,0,cos0) = — 



i(l — rfJo) + ^(3roQ — 1) cos'^ 6 — %/2Re{rjQ} sin 26 cos — rjli sin^ 6 cos 20 



- €Cos2$[rl-^ sin^ 61 + t-qo cos^ 6 - \/2Re{rJo} sin 261 cos - rl_-^ sin^ 6* cos 2</)j 

— e sin 2<?^\/2Ini{riQ} sin 26^ sin (j) + Im{ri_i} sin^ 6^ sin 2(f>^ 

+ v/2e(l + e) cos^(^r-ii sin^ 9 + roo cos^ 5 - \/2Re{rio} sin 261 cos - ri_i sin^ 61 cos 2(^^ 
+ \/2e{l + e) sin<?(^\/2Im{rio} sin 26) sin + Im{r?_i} sin^ sin 2(j!)^ 



W^(<2i,(?i,cos6') = — sjl - e2(^/2Im{r?o}sin26lsin(^ + Im{r?_l}sin^6lsin2^i 

+ \/2e{l - e) cos^(^\/2Im{r[o} sin 261 sin 4> + Im{r-[_i} sin^ S sin 2<^ 
+ ^2e(l — e) sin ^ ^rf i sin^ 6 + roo cos^ 6^ — \/2Re{rfo} sin 2^ cos (f) — r\_i sin' 6 cos 2(j!) 



(25) 



(26) 



4.2.1 s-Channel Helicity Conservation 

It was observed [Tl[2] that, at small the amplitudes 
with a helicity flip in the transition 7* -> F (Ay 7^ A..^) 
are much smaller than those for diagonal transitions where 
the helicity of the vector meson is equal to that of the 
virtual photon (Ay = A^). This behavior is controlled by 
the small parameter 



(27) 



M 



An expansion in /3 corresponds to the twist decomposition. 
A theoretical analysis in the framework of pQCD was per- 
formed [501131] for the amplitudes of the process in Eq. ([T]). 
It was shown [T^ that only the dominant amplitude Too 
receives twist-2 contribution, while all other amplitudes 
contain only higher-twist contributions. The amplitudes 
2^11, Tqi, Ti_i are suppressed by /? while the amplitude 
Tio is suppressed by 0^ . SCHC and the twist decomposi- 
tion lead to the following hierarchy of NPE amplitudes at 
small |<'| and large Q^: 



where M is a typical hadronic mass, usually taken to be 
the nucleon mass. This dominance of the 7* — ^ tran- 
sition that is diagonal with respect to helicity is called 
the s-channel helicity conservation (SCHC) approxima- 
tion. Furthermore, it was shown [40] that for small \t'\ 
and > 2M^, NPE amplitudes with nucleon helicity flip 
are suppressed compared to amplitudes describing diago- 
nal transitions with A at = A^. The same small parameter 
given by Eq. ([?7]) controls this suppression. 

As only an unpolarized target is considered here, there 
is no linear contribution of nucleon- helicity-flip amplitudes 
to the relevant SDMEs [31[S]. The fractional contribution 
of NPE amplitudes with nucleon helicity flip is suppressed 
by a factor a^, so that we need to consider only NPE 
amplitudes with Ajv = A^ if we neglect terms of order . 

The SCHC NPE amplitudes arc Too and Tn. The hehcity- 
flip amplitudes Tqi and Tio are proportional to the small 
factor a, while the double-helicity-flip amplitude Ti_i is 
proportional to . 



4.2.2 Twist decomposition 

At asymptotically large photon virtuality, all amplitudes 
can be decomposed into power series of another small pa- 
rameter 



Mv 



iTooP » |Tn|2 » iToip » iTioP - |Ti_i|2 



(29) 



Such a hierarchy was observed in the HERA collider ex- 
periments [Tl[2ni[Ml[21]- 



4.2.3 Asymptotic Behavior of Amplitudes at High Energy 

At high energy, a third small parameter can be defined: 
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M 
W 



(30) 



(28) 



The asymptotic behavior of amplitudes at large W and 
small |t| was studied experimentally, providing informa- 
tion that led to the development of Regge phenomenol- 
ogy |351I36| . The fractional contribution of various ampli- 
tudes to SDMEs can be estimated by applying the formula 
for the amplitude of the exchange of a reggeon R at small 
1^1 and large W, F (x (VF^/M^)"^^*) , where asit) is the 
Regge trajectory (the spin of the exchanged reggeon). As 
a rough estimate, we assume that the pomeron intercept 
ap(0) ~ 1. Hence the amplitude of pomeron exchange is 
proportional to the factor W'^/M'^. Secondary reggeons 
with natural parity have ajj(O) sa 0.5 which results in 
F (X W/M. Intercepts for reggeons with unnatural par- 
ity are negative in accordance with the results of a Regge 
phenomenology analysis |35ll36j of experimental data on 
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soft scattering of hadrons. Therefore, the Regge factor is 
less than (VF^/M^)° = 1. In QCD, exchanges of secondary 
reggeons and reggeons with unnatural parity correspond 
to quark- ant iquark exchanges, while pomeron exchange 
corresponds to two-gluon exchange. In other words, the 
ratio of amplitudes of quark-antiquark exchanges with 
natural parity to the amplitude of pomeron exchange is 
proportional to the parameter 7. while the ratio of ampli- 
tudes of quark-antiquark exchanges with unnatural parity 
to the amplitude of pomeron exchange is proportional to 
7^ . Since there is no interference between contributions of 
NPE and UPE amplitudes to SDMEs measured on unpo- 
larized targets, the fractional contribution of UPE ampli- 
tudes is suppressed by a factor 7^. This explains why only 
NPE amplitudes survive at HERA collider energies. 

In the context of singie-reggeon exchange, the isospin 
1=1 reggeon (p, a2, ...) contribution to the amplitude of 
vector-meson production on the proton is of opposite sign 
to that on the neutron. In contrast, the / = reggeon 
(pomeron, w, /2, ...) contribution to the amplitude is of 
the same sign for both proton and neutron. Hence, the 
results on SDMEs and amplitude ratios should be differ- 
ent for proton and deuteron if / = 1 and 1 = contri- 
butions to the amplitude interfere. By comparing proton 
and deuteron results, the fractional contribution of / = 1 
reggeons to the extracted amplitude ratios can be esti- 
mated. 



4.2.4 Hierarchy of Amplitudes in the Kinematic Region of 
HERMES 

In the HERMES kinematic region, the parameter /3 is 
larger than 0.3 and the relative sizes of SDMEs measured 
in exclusive production can be explained by the follow- 
ing hierarchy [TSj : 

\Too\' - iTiip » iC/iip > iToip » iTioP 

^|ri_ip>|[/oi|2^|C/ioP~|C/i-iP, (31) 

with the two largest amplitudes being Too and Tu. The 
abbreviated notation U\^x^, where 

Px.xj' ^ |C/A,iA._i|' + |C^A.iA,i|' (32) 

was introduced because, for UPE amplitudes, it is im- 
possible to prove the dominance of those without nucleon 
helicity flip over those with helicity flip, in contrast to 
the NPE amplitudes. As shown in Eq. (PT|) . the moduli 
of all the UPE amplitudes except Uu are smaller than 
those of the NPE amplitudes. The numerical estimate 
|t^ii|V(|Tii|2-He|TooP) ~ 0.11 ±0.04 obtained in Ref. [HI 
shows that the modulus of \Uii\ is even larger than that 
of |Toi| (and the result of the fit done in the present work 
confirms this result) and the contribution of \Uii\ to the 
SDMEs cannot be neglected. This contribution is sup- 
pressed by a factor 7^ which is smaller than 0.01 in the 
HERMES kinematic region 3.0 GeV < < 6.5 GeV. If 
the UPE amplitude Uu is due to pion exchange, its con- 
tribution may be significant in the HERMES kinematic 



region because of the large pion-nucleon coupling constant 

Contributions from small amplitudes can be apprecia- 
ble if they are multiplied by the largest amplitudes Too or 
Til in the numerators of SDME formulas. As the small am- 
plitudes Toi, Tio, and Ti-i are multiplied by the largest 
amplitudes (sec Ref. [3] and Appendix A of Ref. [18] ). 
five complex NPE amplitudes have to be considered in 
total, i.e. the four ratios Tu/Tqq, Tqi/Tqq, Tiq/Too, and 
Ti_i/Too. Concerning the ratio C/u/Too, only |i7ii|/|Too| 
can be determined as there is no interference between UPE 
and NPE amplitudes for the case of an unpolarized tar- 
get. The contribution of all other UPE amplitudes can 
be neglected [T^ . In total, there are nine independent free 
parameters to be determined when fitting the angular dis- 
tribution using the amplitude method. 

4.3 Basic Formulas of the Amplitude Method 

The exact formulas for SDMEs expressed in terms of he- 
licity amplitudes are presented in Appendix A of Ref. [TB] . 
In order to rewrite SDMEs in terms of the ratios of sig- 
nificant helicity amplitudes, we neglect the contributions 
from NPE nucleon helicity-flip amplitudes (using Eq. (ITSl) ') 
and from all UPE amplitudes except Uu. Then, dividing 
both the numerators and the denominator in the exact 
formulas for SDMEs by |TooP, we arrive at approximate 
expressions for SDMEs in terms of certain amplitude ra- 
tios: 



'00 — 


{e+\toi\^}/N, 




(33) 


Re{r?on ^ 


Re{eiio + ^toi(^ 


ii-ti-i)*}/VV, 


(34) 


^1-1 — 
"11 — 
'00 — 


Rc{-e|iio|'+ti 
Re{ii_iiti}/iV, 
-\toi\VN, 




(35) 
(36) 
(37) 


Re{r}o} ^ 


iRc{-ioi(tii - 


h-in/N, 


(38) 


^1-1 — 






(39) 


Im{r^o} ~ 






(40) 


Im{^i-i} " 


l{-\hi\' + \h- 


i\' + \uii\'}/N, 


(41) 


^5 _ 
'11 — 


-^Re{iio(tii - 




(42) 


'00 — 


V2Re{toi}/N, 




(43) 


Rc{r%} ^ 


-^Re{2tiotSi + 


(ill -ii-i)}/iV, 


(44) 


^1-1 — 


--^Re{-tio(tii 




(45) 


lm{r%} ^ 


-^Rc{tii + ti 
v8 




(46) 


Im{rti} ^ 


--^Re{iio(iii + 




(47) 
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Ini{r3j ~ --Im{ioi(iii +ii-i)*}/A^, (48) 

Ini{rti} ^ -Im{ti_itti}/7V , (49) 

lm{rlo} - ^M^ii + , (50) 

Im{rl_i} ~ -Llm{tio(<ii + <i-i)*}/iV, (51) 



' 00 



1 



r«i~ ^Im{tio(tii -ti-i)*}/iV, 



V2lm{toi}/N, 



Rc{r%} ~ ^Im{-2iio<Si (54) 



1 
1 

75^ 



rr_i 0^ -^Imitwihi - ti-i)*}/iV 



(52) 
(53) 



(55) 



where the normahzation factor TV = A^/|TooP is defined 
in accordance with Eqs. (|19I) - (PT|) by 



neghgible, while the rapid (60- 180s) reversal of the polar- 
ization direction minimized polarization bias due to de- 
tector effects. 

HERMES was a forward spectrometer [5T] in which 
both scattered lepton and produced hadrons were detected 
within an angular acceptance of ±170 mrad horizontally, 
and ±(40— 140) mrad vertically. The tracking system had 
a momentum resolution of about 1.5% and an angular res- 
olution of about 1 mrad. Lepton identification was accom- 
plished using a transition-radiation detector, a preshower 
scintillator counter, and the electromagnetic calorimeter. 
The particle identification system included a gas threshold 
Cerenkov counter, which was replaced in 1998 by a dual- 
radiator ring-imaging Cerenkov detector |45| . Combining 
the responses of these detectors in a likelihood method led 
to an average lepton identification efficiency of 98% with a 
hadron contamination of less than 1%. The event sample 
used in this analysis is exactly the same as that used in 
Ref. UHl. 



Nt^ |iii|' + |toi|' + |ii-i|' + |wii|' 



Nl ~ l + 2|tio|2. 
Here txyX^ = Tx^x^/Toa and |mii| 



(56) 

|C/ii|V|7^oo|' with 
|C/ii|^ defined in Eq. (|32|) . There are nine independent 
real functions in Eqs. p3ll56|) . namely: Re(tii), Im(tii), 
Re(toi), Im(<oi), Re(tio), Im(iio), Re(<i_i), Im(ti_i), and 



Equations p3ll56p are the basis for the extraction of 
the helicity amplitude ratios from the measured angular 
distributions and also for the calculation of SDMEs from 
amplitude ratios (see Sec. 16. 5p . 



6 Extraction of Amplitude Ratios 

The measurement of the angular distribution of the scat- 
tered electrons/positrons and the pions from the decay 
— > 7r+7r~ is the basis for the extraction of spin density 



matrix elements 



in the SDME method and of he- 



Xv 

licity amplitude ratios ^Av-a^ and |wii| in the amplitude 
method. For a polarized lepton beam and an unpolarized 
target, the angular distribution of the pions is given by 
relations ([Mll^ in Sec. (see also Refs. 011]), where 
the SDMEs are expressed in terms of helicity amplitude 
ratios according to Eqs. P3II56|) . 



5 The HERMES Experiment 

A detailed description of the HERMES experiment at 
DESY can be found in Ref. [H] . The data analyzed in this 
paper were collected between the years 1996 and 2005. 
A longitudinally polarized positron or electron beam of 
27.6 GeV was scattered from pure hydrogen or deuterium 
gas targets internal to the HERA lepton storage ring. The 
lepton beam was transversely polarized due to an asym- 
metry in the emission of synchrotron-radiation [32] in the 
field of the dipolc magnets. The transverse beam polariza- 
tion was transformed locally into longitudinal polarization 
at the interaction point by spin rotators located upstream 
and downstream of the HERMES apparatus. The helic- 
ity of the beam was typically reversed every two months. 
For both positive and negative beam helicities, the beam 
polarization was continuously measured by two Compton 
polarimeters (151144] . The average beam polarization for 
the hydrogen (deuterium) data set was 0.45 (0.47) after 
requiring 0.15 < Pb < 0.8 in the analysis process, and 
the fractional uncertainty of the beam polarization was 
3.4% (2.0%) ^MM. Some of the data were collected with 
longitudinally or transversely polarized targets. However, 
the time-averaged polarization of the polarized targets was 



6.1 Binned Maximum Likeliliood Method 

Amplitude ratios are extracted with the same binned max- 
imum likelihood method as in the previous SDME analysis 
(see Sec. 6 of Ref. [18] ). In brief, the amplitude ratios in 
each of the kinematic bins are obtained by minimizing the 
difference between the 3-dimensional (cos6', </>, ^) angular 
distribution of the experimental events and that of a sam- 
ple of fully reconstructed Monte Carlo events. The latter 
are generated isotropically in (cos 9, </>, <?) using the rhoMC 
generator [TflH5] for exclusive p° simulated in the instru- 
mental context of the spectrometer, and passed through 
the same reconstruction chain as the experimental data. 
The generated Monte Carlo events are iteratively reweigh- 
ted with the angular distribution given by Eqs. (j24ll26|) . 
with the amplitude ratios treated as free parameters. 

The minimization itself and the uncertainty calcula- 
tion are performed using the MINUIT package [IS]- More 
details can be found in Ref. [T8] . 

6.2 Results on Amplitude Ratios in Bins of and —t' 

Ratios of helicity amplitudes are extracted in a two-dimen- 
sional {Q^, —t') binning from the same HERMES proton 
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bin 


(0'>, GeV^ 


{-*'), GeV^ 


bin 




{-*'), GeV^ 


qiti 


0.817 


0.019 


93*1 


1.658 


0.019 


91*2 


0.823 


0.068 


93*2 


1.660 


0.068 


qiti 


0.821 


0.146 


93*3 


1.663 


0.146 


qiti 


0.815 


0.280 


93*4 


1.663 


0.282 


(72*1 


1.184 


0.019 


94*1 


2.996 


0.019 


92*2 


1.188 


0.068 


94*2 


3.056 


0.068 


(72*3 


1.189 


0.145 


94*3 


3.076 


0.146 


92*4 


1.188 


0.282 


94*4 


3.134 


0.284 



Table 1. Mean values of kinematic variables for 16 bins. The limits of the gi, 92, 93, and 94 bins for are 0.5; 1.0; 1.4; 2.0; 
7.0 GeV^ while those of the *i, *2, *3, and *4 bins for -*' are the following: 0.0; 0.04; 0.10; 0.20; 0.40 GeV^ 



and deuteron data sets as were used for the SDME analysis 
of Ref. [18]. The four bins are denoted as qi, q2, 93, and 
94 defined by the limits 0.5, 1.0, 1.4, 2.0, 7.0 GeV^. The 
four —t' bins are denoted as ii, *2, ^3, and *4 defined by the 
limits 0.0, 0.04, 0.10, 0.20, 0.40 GeV^. The mean values 
of and — *' in each bin are presented in Tab. [TJ The 
results of the extraction of the amplitude ratios in 16 bins 
are presented in Tab. [2] for the proton data and in Tab. [3] 
for the deuteron data. In every {Q^, —t') bin, the nine free 
parameters are obtained from a fit to the 3-dimensional 
angular distribution without subtracting the background 
from semi- inclusive deep- inelastic scattering (SIDIS). The 
SIDIS background under the exclusive peak, i.e., in the 
region —1 GcV < AE < 0.6 GeV, is estimated using the 
PYTHIA Monte Carlo generator [SD] (see Ref. [H] for 
more details) and its effect is assigned as a systematic 
uncertainty. 



6.3 Testing the Extraction Method 

The self-consistency of the extraction of amplitude ra- 
tios is tested with rhoMC Monte Carlo data. The fit re- 
sults obtained for the nine free parameters in several {Q^, 
—t') bins spanning the experimental kinematic range are 
used to calculate the 23 SDMEs according to Eqs. 
[5^ . These SDMEs are then used as input for a Monte 
Carlo simulation of the angular distribution described in 
Eqs. (j24ll26|) . This angular distribution is used to extract 
the nine free parameters again. The resulting amplitude 
ratios are found to be consistent with the input amplitude 
ratios within statistical uncertainties. 



uncertainty due to the background is chosen to be the dif- 
ference between the amplitude ratios of these two sets. The 
other main uncertainty arises from the reliance on Monte 
Carlo simulations to perform the extraction of amplitude 
ratios from data. There are uncertainties in the parame- 
ters of the description of production in the Monte Carlo 
that propagate through to the amplitude-ratio values pre- 
sented in this paper. The parameterization of the total 
electroproduction cross section in rhoMC is chosen in the 
context of a vector-meson-dominance model [T] that incor- 
porates a propagator-type dependence, and also con- 
tains a dependence of the ratio of the longitudinally- 
polarized virtual photon and the transversely-polarized 
virtual photon cross-sections. As the HERMES spectrom- 
eter acceptance depends on Q^, different input parameters 
result in slightly different reconstructed isotropic angular 
distributions. The corresponding systematic uncertainty 
of the resulting helicity amplitude ratios is obtained by 
varying these parameters within one standard deviation in 
the total uncertainty of the parameters given in Refs. p71 
H51H71H5] . The effect of a possible kinematic dependence 
of the *' slope of the cross section in the rhoMC generator 
is found to be negligible in the kinematic region of our 
experiment. 

A third main contribution to the systematic uncer- 
tainty arises in the amplitude method as a result of the 
neglect of the small NPE nucleon helicity-flip amplitudes 
and all the UPE amplitudes except U^^^ij^i. As shown 
in Appendix [21 the real and imaginary parts of the de- 
viation (5*Av-At and also of the obtained amplitude 
ratios from the true values tx^x^ and |wii|, respectively, 
can be estimated using the relations 



6.4 Systematic Uncertainties 

In the extraction of SDMEs in Ref. [15], two main con- 
tributions to the total systematic uncertainty of the re- 
sults were discussed. One was the uncertainty in the back- 
ground contribution to the signal in the region of the ex- 
clusive peak in the AE distribution. The background was 
simulated with the PYTHIA Monte Carlo generator. In 
the present analysis, the amplitude ratios are determined 
with and without background subtraction. The systematic 



\5Re{tx,.x-,)\^^\tx.x^, (57) 

|5Im(iA.Aj| = ^ItA.Aj, (58) 

5\uii\ = ^\uii\. (59) 

These equations are used to calculate the systematic un- 
certainty for the amplitude method described in this pa- 
per. 
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Ratio 


qiti 


9l*2 


91*3 


91*4 


Re(tii) 


0.975 ±0.121 ±0.297 


0.961 ±0.101 ±0.213 


1.363 ±0.140 ±0.309 


1.037 ±0.176 ±0.230 


Im(iii) 


0.542 ±0.187 ±0.052 


0.285 ±0.145 ±0.174 


0.082 ±0.285 ±0.104 


0.783 ±0.147 ±0.213 


Re(toi) 


0.025 ± 0.044 ± 0.056 


0.112 ±0.036 ±0.047 


0.214 ± 0.064 ± 0.050 


0.182 ±0.048 ±0.116 


Im(toi) 


-0.098 ±0170 ±0.211 


0.111 ±0.111 ±0.136 


0.315 ±0.120 ±0.073 


0.107 ±0.108 ±0.107 


Re(tio) 


0.037 ±0.038 ±0.050 


-0.049 ± 0.039 ± 0.024 


0.024 ± 0.039 ± 0.026 


-0.010 ± 0.056 ±0.043 


Im(tio) 


-0.067 ±0.064 ±0.156 


-0.066 ± 0.068 ± 0.009 


0.005 ± 0.069 ± 0.023 


0.050 ±0.067 ±0.021 


Re(ti-i) 


-0.110 ±0.042 ±0.045 


-0.021 ±0.037 ±0.079 


-0.037 ± 0.043 ±0.043 


0.020 ± 0.056 ± 0.040 


Im(ti-i) 


0.178 ±0.087 ±0.349 


-0.147 ±0.092 ±0.055 


-0.124 ±0.100 ±0.079 


-0.172 ± 0.079 ±0.104 


|«ii| 


0.329 ±0.070 ±0.021 


0.424 ± 0.049 ± 0.032 


0.391 ±0.063 ±0.124 


0.357 ±0.077 ±0.070 


Ratio 


92*1 


92*2 


92*3 


92*4 


Re(tii) 


1.138 ±0.143 ±0.021 


0.996 ±0.106 ±0.084 


1.079 ±0.094 ±0.135 


0.925 ±0.097 ±0.139 


Im(tii) 


0.282 ±0.268 ±0.212 


0.386 ±0.125 ±0.039 


0.342 ±0.162 ±0.159 


0.289 ±0.156 ±0.141 


Re(toi) 


0.044 ± 0.051 ± 0.009 


0.116 ±0.037 ±0.061 


0.113 ±0.035 ±0.030 


0.250 ± 0.043 ± 0.075 


Im(toi) 


0.073 ±0.146 ±0.137 


0.327 ± 0.099 ± 0.083 


-0.009 ±0.141 ±0.069 


0.265 ±0.107 ±0.055 


Re(tio) 


0.001 ±0.054 ±0.011 


-0.054 ±0.034 ±0.033 


0.032 ± 0.030 ± 0.015 


-0.065 ± 0.032 ±0.020 


Im(tio) 


-0.081 ±0.073 ±0.154 


0.064 ± 0.073 ± 0.058 


0.006 ± 0.070 ± 0.027 


0.063 ± 0.073 ± 0.057 


Re(ti_i) 


0.001 ±0.045 ±0.015 


0.023 ±0.040 ±0.018 


-0.017 ± 0.036 ±0.021 


-0.073 ± 0.035 ±0.019 


Im(ti_i) 


-0.041 ±0.137 ±0.271 


0.002 ±0.078 ±0.019 


-0.084 ±0.107 ±0.024 


-0.099 ± 0.069 ±0.034 


'till 


0.499 ± 0.040 ± 0.039 


0.429 ±0.042 ±0.081 


0.418 ±0.046 ±0.119 


0.359 ± 0.055 ± 0.055 


Ratio 


93*1 


93*2 


93*3 


93*4 


Re(tii) 


1.029 ±0.104 ±0.039 


1.088 ± 0.098 ± 0.083 


0.878 ± 0.079 ± 0.202 


0.975 ± 0.096 ± 0.228 


Im(tii) 


0.257 ±0.131 ±0.151 


0.479 ±0.111 ±0.050 


0.513 ±0.102 ±0.385 


0.280 ±0.176 ±0.365 


Re(toi) 


-0.052 ±0.039 ±0.052 


0.008 ± 0.032 ± 0.060 


0.113 ±0.035 ±0.081 


0.194 ±0.044 ±0.064 


Im(toi) 


0.357 ±0.127 ±0.480 


0.100 ±0.100 ±0.004 


0.123 ±0.092 ±0.275 


-0.100 ±0.108 ±0.406 


Re(tioO 


0.074 ± 0.028 ± 0.021 


0.043 ±0.028 ±0.011 


-0.007 ±0.034 ±0.032 


-0.032 ± 0.038 ± 0.036 


Im(tio) 


-0.051 ±0.067 ±0.042 


-0.129 ±0.058 ±0.031 


0.075 ±0.053 ±0.173 


0.025 ±0.113 ±0.120 


Re(ti_i) 


-0.013 ±0.038 ±0.021 


-0.019 ±0.035 ±0.037 


0.034 ± 0.037 ±0.101 


-0.044 ± 0.036 ±0.065 


Im(ti_i) 


0.108 ±0.090 ±0.236 


0.100 ±0.071 ±0.031 


-0.186 ±0.071 ±0.071 


0.037 ±0.098 ±0.185 


|«ii| 


0.423 ± 0.055 ± 0.048 


0.323 ± 0.068 ± 0.084 


0.346 ± 0.056 ± 0.085 


0.445 ±0.050 ±0.191 


Ratio 


94*1 


94*2 


94*3 


94*4 


Re(tii) 


0.723 ±0.071 ±0.053 


0.706 ± 0.068 ± 0.039 


0.582 ±0.074 ±0.156 


0.650 ±0.063 ±0.117 


Im(fii) 


0.570 ±0.071 ±0.025 


0.657 ±0.061 ±0.066 


0.583 ±0.110 ±0.120 


0.488 ±0.067 ±0.119 


Re(toi) 


0.062 ± 0.031 ± 0.007 


0.121 ±0.031 ±0.006 


0.190 ±0.054 ±0.021 


0.282 ± 0.036 ± 0.060 


Im(toi) 


0.011 ±0.069 ±0.019 


0.084 ±0.079 ±0.019 


-0.129 ±0.0.57 ±0.154 


-0.099 ±0.081 ±0.044 


Re(tio) 


-0.013 ±0.030 ±0.011 


-0.046 ± 0.035 ± 0.008 


0.095 ±0.025 ±0.118 


0.007 ±0.031 ±0.012 


Im(iio) 


0.010 ±0.052 ±0.010 


0.050 ±0.045 ±0.015 


-0.181 ± 0.060 ±0.258 


0.003 ± 0.044 ± 0.058 


Re(ti-i) 


0.000 ± 0.034 ± 0.007 


0.065 ± 0.034 ± 0.007 


0.027 ± 0.040 ± 0.033 


0.024 ±0.031 ±0.011 


Im(ti_i) 


0.029 ± 0.056 ± 0.005 


-0.068 ±0.048 ±0.015 


-0.038 ± 0.040 ±0.142 


-0.036 ±0.051 ±0.066 


|mii1 


0.451 ± 0.049 ± 0.024 


0.306 ± 0.061 ± 0.032 


0.383 ±0.067 ±0.097 


0.380 ± 0.044 ± 0.091 



Table 2. Ratios of lielicity amplitudes in and —t' bins for proton data. The mean values of and —t' for all 16 bins 
are given in Tab.[T] The notations tkj and |?iii| are used for the amplitude ratios Tkj /Too and jf/ii/Too], respectively. The first 
uncertainties are statistical, the second are systematic. 



The three main systematic uncertainties are added in 
quadrature to give the total systematic uncertainty pre- 
sented in Tabs. and [31 

The effect of radiative corrections onto the values of 
SDMEs extracted was shown to be negligible [15] . This 
is mainly due to the exclusivity cut AE < 0.6 GeV that 



excludes photons radiated with more than 0.6 GeV. The 
amplitude ratios in this work are extracted using exactly 
the same data sets and cuts, by a very similar fit com- 
paring the shapes of the 3-dimensional experimental and 
simulated angular distributions. Hence we conclude that 
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Ratio 


ijiti 


91*2 


9l*3 


9l*4 


Re(tii) 


0.860 ±0.078 ±0.163 


1.256 ± 0.095 ± 0.055 


1.220 ±0.113 ±0.145 


1.197 ±0.124 ±0.218 


Im(tii) 


0.509 ±0.118 ±0.203 


0.304 ±0.159 ±0.049 


0.192 ±0.228 ±0.114 


0.222 ± 0.320 ± 0.208 


Re(toi) 


-0.011 ± 0.029 ±0.027 


0.020 ± 0.032 ± 0.069 


0.174 ±0.044 ±0.080 


0.151 ±0.052 ±0.042 


Im(toi) 


0.023 ± 0.076 ± 0.070 


0.088 ± 0.076 ± 0.035 


0.239 ±0.112 ±0.086 


0.222 ±0.143 ±0.164 


Re(tio) 


-0.011 ± 0.029 ±0.020 


0.002 ±0.027 ±0.014 


0.025 ± 0.032 ± 0.048 


0.017 ±0.033 ±0.019 


Im(tio) 


0.020 ± 0.049 ± 0.044 


-0.017 ±0.063 ±0.047 


-0.048 ±0.083 ±0.135 


0.014 ±0.076 ±0.095 


Re(ti-i) 


-0.017 ± 0.029 ±0.018 


0.056 ± 0.030 ± 0.051 


-0.052 ± 0.032 ± 0.039 


0.027 ± 0.050 ± 0.078 


Im(ti_i) 


-0.062 ± 0.064 ±0.020 


-0.098 ± 0.078 ± 0.039 


-0.093 ±0.093 ±0.153 


-0.228 ±0.105 ±0.054 




0.364 ±0.037 ±0.090 


0.343 ± 0.046 ± 0.037 


0.419 ± 0.048 ± 0.020 


0.388 ±0.060 ±0.105 


Ratio 


92*1 


92*2 


92*3 


92*4 


Re(tii) 


0.984 ±0.087 ±0.048 


0.888 ±0.080 ±0.126 


1.119 ±0.074 ±0.129 


0.990 ±0.063 ±0.140 


Im(tii) 


0.341 ±0.106 ±0.054 


0.778 ±0.078 ±0.153 


0.252 ±0.105 ±0.163 


0.039 ±0.142 ±0.482 


Re(toi) 


0.042 ± 0.023 ± 0.022 


0.071 ±0.028 ±0.028 


0.109 ±0.025 ±0.037 


0.262 ± 0.031 ± 0.043 


Im(toi) 


0.289 ±0.130 ±0.084 


-0.062 ±0.069 ±0.158 


0.136 ±0.074 ±0.035 


0.070 ± 0.086 ± 0.064 


Re(tio) 


-0.006 ± 0.027 ±0.020 


0.039 ± 0.025 ± 0.038 


0.039 ±0.019 ±0.012 


0.022 ± 0.029 ± 0.048 


Im(tio) 


0.029 ±0.101 ±0.052 


-0.009 ± 0.034 ± 0.013 


-0.087 ± 0.044 ± 0.007 


-0.201 ±0.065 ±0.171 


Re(ti-i) 


0.013 ±0.046 ±0.017 


0.021 ± 0.028 ± 0.060 


-0.077 ±0.007 ±0.012 


-0.041 ± 0.028 ± 0.039 


Im9ti_i) 


0.063 ±0.100 ±0.034 


-0.042 ± 0.046 ± 0.042 


0.030 ± 0.070 ± 0.012 


-0.071 ±0.072 ±0.106 


Mil 


0.395 ± 0.035 ± 0.066 


0.262 ± 0.050 ± 0.048 


0.401 ±0.036 ±0.134 


0.373 ± 0.048 ± 0.085 


Ratio 


93*1 


93*2 


93*3 


93*4 


Re(tii) 


0.787 ± 0.059 ± 0.052 


0.819 ±0.065 ±0.067 


0.839 ± 0.076 ± 0.076 


0.835 ±0.060 ±0.128 


Im(tii) 


0.487 ± 0.072 ± 0.081 


0.518 ±0.085 ±0.079 


0.553 ±0.088 ±0.111 


0.203 ±0.097 ±0.129 


Re(toi) 


0.028 ± 0.024 ± 0.030 


0.094 ± 0.025 ± 0.053 


0.129 ±0.028 ±0.017 


0.222 ± 0.028 ± 0.034 


Im(toi) 


0.102 ±0.080 ±0.021 


0.001 ±0.091 ±0.042 


-0.017 ±0.078 ±0.072 


0.080 ± 0.072 ± 0.040 


Re(tio) 


-0.006 ± 0.024 ±0.003 


-0.007 ± 0.023 ± 0.020 


0.027 ±0.028 ±0.010 


0.053 ± 0.026 ± 0.038 


Im(tio) 


0.008 ± 0.048 ± 0.018 


-0.025 ± 0.045 ± 0.042 


-0.092 ±0.046 ±0.014 


-0.228 ± 0.058 ± 0.035 


Re(ti-i) 


0.031 ±0.026 ±0.026 


-0.045 ±0.027 ±0.012 


-0.018 ±0.032 ±0.052 


-0.008 ±0.027 ±0.025 


Im(ii_i) 


0.000 ± 0.052 ± 0.037 


0.102 ±0.061 ±0.015 


-0.038 ± 0.065 ± 0.048 


-0.068 ± 0.067 ± 0.047 


kill 


0.386 ± 0.034 ± 0.020 


0.402 ± 0.039 ± 0.076 


0.355 ± 0.045 ± 0.023 


0.353 ± 0.050 ± 0.072 


Ratio 


94*1 


94*2 


94*3 


94*4 


Re(tii) 


0.655 ± 0.054 ± 0.044 


0.842 ± 0.068 ± 0.043 


0.768 ± 0.060 ± 0.079 


0.820 ±0.075 ±0.315 


Im(tii) 


0.629 ±0.049 ±0.011 


0.614 ± 0.066 ± 0.052 


0.592 ± 0.059 ± 0.089 


0.614 ± 0.070 ± 0.273 


Re(toi) 


0.060 ± 0.025 ± 0.003 


0.147 ±0.030 ±0.010 


0.144 ±0.029 ±0.018 


0.243 ±0.031 ±0.105 


Im(ioi) 


-0.091 ± 0.075 ±0.011 


-0.065 ±0.067 ±0.018 


-0.173 ±0.076 ±0.028 


0.042 ± 0.078 ± 0.382 


Re(tio) 


0.025 ±0.027 ±0.009 


0.042 ±0.030 ±0.012 


0.073 ±0.031 ±0.016 


-0.014 ±0.034 ±0.063 


Im(tio) 


-0.088 ± 0.038 ±0.014 


-0.148 ±0.047 ±0.010 


-0.131 ±0.041 ±0.017 


-0.025 ± 0.049 ± 0.043 


Re(ti-i) 


-0.009 ± 0.030 ±0.001 


-0.004 ±0.032 ±0.017 


-0.010 ±0.030 ±0.008 


-0.009 ± 0.030 ± 0.046 


Im(ti_i) 


0.016 ±0.038 ±0.006 


-0.024 ±0.051 ±0.026 


-0.050 ± 0.042 ± 0.075 


-0.055 ±0.048 ±0.109 


kill 


0.436 ± 0.038 ± 0.042 


0.422 ± 0.044 ± 0.045 


0.389 ± 0.042 ± 0.070 


0.420 ± 0.045 ± 0.076 



Table 3. Ratios of helicity amplitudes in and —t' bins for deuteron data. The mean values of and —t' for all 16 bins 
are given in Tab.[T] The notations tkj and |?iii| are used for the amplitude ratios Tkj /Too and jf/ii/Too], respectively. The first 
uncertainties are statistical, the second are systematic. 



any possible systematic uncertainty of amplitude ratios 
due to radiative corrections can be safely neglected. 



6.5 Calculation of SDMEs from Extracted Helicity 
Amplitude Ratios 

In the following, the SOME results obtained from SOME 
and amplitude methods are compared. The SDMEs are 
calculated for every (Q^, — *') bin by applying formu- 
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2. Comparison of SDME and amplitude methods. Red squares show the resuhs of the SDME analysis |18) . Blue circles 



(amplitude method) are obtained in the present work from amplitude ratios fitted directly to the three-dimensional angular 
distribution in every Q2 and -t' bin. The proton data at (Q^) = 3 GeV^ {-t') = 0.019 GeV^ are presented. Yellow bands mark 
those SDMEs that are non-zero in the SCHC approximation. Total uncertainties are depicted. 



las p3ll56p and using the values of the real and imagi- 
nary parts of the amplitude ratios given in Tabs. and 
[31 The results of both methods are found to agree over 
the full kinematic range of the experiment. The compar- 
ison is shown for two representative bins, for the proton 
data using the bin 94*1 ((Q^) = 3 GeV^ {-t') = 0.019 
GeV^) in Fig. [5] and for deutcron data using the bin 52^3 
((Q2) = 1.19 GeV2, {-t') = 0.145 GeV^) in Fig. H The 
bands in the figures mark those SDMEs that may be non- 



zero in the SCHC approximation, while the line at zero is 
drawn for those SDMEs that vanish at t' = 0. 

For the unpolarized SDMEs, the total uncertainties 
obtained with the amplitude method are in most cases 
comparable to those obtained with the SDME method. 
An exception seen in Fig. [2] is Tqq which is proportional to 
1^01 P (see Eq. ([57)) ) and hence extremely small at small 
\t'\. For polarized SDMEs, the uncertainties obtained us- 
ing the SDME method are generally much larger than 
those obtained using the amplitude method. This differ- 
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Fig. 3. Comparison of SDME and amplitude methods. Red squares show the results of the SDME analysis [18]. Blue circles 
(amplitude method) are obtained in the present work from amplitude ratios fitted directly to the three-dimensional angular 
distribution in every and —t' bin. The deuteron data at {Q'^) = 1.19 GeV^, {—t') = 0.145 GeV'^ are presented. Yellow bands 
mark those SDMEs that are non-zero in the SCHC approximation. Total uncertainties are depicted. 



ence reflects a major difference in the methods themselves. 
In the SDME method, polarized and unpolarized SDMEs 
are fitted as independent free parameters. In this case, the 
error bars of the polarized SDMEs are inflated due to the 
factors Pb = 0.47 and ^/^^ w 0.45 (see Eqs. IQMB)- 
However, Eqs. p3ll56p show that the polarized SDMEs 
depend on amplitude ratios constrained by both polar- 
ized and unpolarized data. Hence the polarized SDMEs 
calculated from the amplitude ratios have uncertainties 
comparable to those of the unpolarized SDMEs. Never- 



theless, the polarized data serves the important function 
of constraining the sign of the imaginary parts of ampli- 
tude ratios. 
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7 Kinematic Dependences of Amplitude 
Ratios 

7.1 Predictions for Asymptotic Behavior 

At small |t'|, the behavior of NPE amplitudes T\y\_^ can 
be described [5] by 



the ratios tn, ioi- a-nd ti_i (tio) are odd (even) functions 
in Q. The fit functions will be chosen in agreement with 
this property of amplitude ratios. Whenever the chosen 
fit function does not agree with a pQCD prediction, the 
resulting curve is shown not by a solid but a dash-dotted 
line. 



(60) 7.2 Kinematic Dependence of Tu/Ti 



which refiects angular momentum conservation in the pro- 
cess in Eq. ([1]) . Equation (pO)) shows that Too behaves as a 
constant at smaU Hence each ratio tAv-A-, T'ai-a^/Tqo 
has the same asymptotic behavior at — > as the ampli- 
tude Ta^^a^. 

The asymptotic behaviour of the amplitudes of vector- 
meson electroproduction at large and small xb was 
considered in Refs. [501151] in the framework of pQCD, 
using the approximation of leading logarithms in which 
the logarithms ln(l/a;s) and In(Q^) are large, and the 
amplitudes are considered to be nearly imaginary. It was 
predicted [SUIISTIIS^ that the amplitude ratios should have 
the following behavior at large and small xb- 



hi = Jii/Joo oc — — , 



Toi/Tc 



00 oc 



iio = Tiq/Tqq oc 



[-1-^1 



= Ti-i/Too oc 



-t'Mv 



Ml 



9i 



(61) 
(62) 
(63) 

(64) 



The functions Ci and C2 in Eq. are dimensionless and 
the parameter /i is between 0.7 GeV and 1.0 GeV. It was 
noted in Refs. [SOll^Tll^^] that Ci is essentially a constant 
while C2 is a ratio whose numerator contains the gluon 
transversity distribution [53] Gt{xb, IJ''^) at the scale 
and the denominator contains the unpolarized gluon den- 
sity G{xb, Qv) rneasured at the conjectured [31] scale Qv 
for vector-meson electroproduction in the non-asymptotic 
region, with 



(65) 

As explained in Sec. 14.2.^1 the twist-decomposition pa- 
rameter j3 = MyjQ is not really small in the kinematic 
region of the HERMES experiment. Therefore, HERMES 
results on amplitude ratios are expected to be not al- 
ways in agreement with the pQCD predictions given by 
Eqs. (|6m64p . In the case of a disagreement, we use a pa- 
rameterization that does not contradict general principles 
and describes the amplitude ratios with reasonable psr 
degree of freedom. Any NPE amplitude ratio has to obey 
Eq. (|60p . On very general principles, at finite Q the am- 
plitudes Til, Tbi, and Ti_i are even functions in Q. In 
contrast, the amplitudes Too and Tio are odd functions 
in Q due to the extra factor l/Q in Eq. (jl2p . Therefore, 



00 



The amplitudes Tqo and Tn describe the diagonal helic- 
ity transitions 7^ — >■ and 7^ — > p\ respectively, and 
are the largest amplitudes of meson production in the 
HERMES kinematic region. The dependence of the ex- 
tracted amplitude ratio in in four —t' bins is presented in 
Fig. m for the proton (left) and the deuteron (right). The 
points correspond to the amplitude ratios extracted from 
the data in 16 bins and given in Tabs. [U and [3] In order 
to test the predictions of Refs. [301131] . the Q"^ dependence 
of the ratio Re(<ii) was fitted in every —t' bin with the 
function 



Re(ii 



Q' 



(66) 



For the four — t' bins, the values of the parameter a are 
found to be consistent with each other within experimen- 
tal accuracy. Therefore, data in all —t' bins are fitted si- 
multaneously to improve the statistical accuracy in deter- 
mining the value of a. The results of the fit for proton 
and deuteron data are presented in Tab. [3] and shown in 
Fig. [U for proton (left) and deuteron (right). The figure 
shows that the dependence of Re(iii) is well described 
by the parameterization given by Eq. ([66]) . 

The dependences of the other amplitude ratios, which 
are described in the following, are examined in a sim- 
ilar manner as was done for Re(tii). In every case, it 
was checked that the respective four datasets containing 4 
points each can be combined into one dataset containing 
16 points. When proton and deuteron results are consis- 
tent, the data are combined and 32 points are used in the 
fit. All the fits to the or t' dependences use 16 data 
points for the proton or deuteron separately and 32 data 
points for the combined data sets. Only for the sake of 
a more clear representation of the kinematic dependence 
of amplitude ratios in one variable (or t'), we average 
the values of the amplitude ratios T = ^AvA.^ or |itii| over 
four —t' (or Q^) bins at the same value of (or t') using 
the standard relations for the mean value (T) and squared 
total uncertainty ((5T)^: 



(r) = 
{8rf = 



E™i/K„)2 ' 

1 

E™i/K„)2 ■ 



(67) 
(68) 



Here, Uj denotes the total uncertainty of in the j bin 
and 8T is the total uncertainty of the averaged ratio T ■ 
The same formula ([05]) is applied for the calculation of the 
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Fig. 4. The dependence of Re(Tii/Too) for proton (left panel) and deuteron (right panel) data, showing the result for the 
16 (Q^, -t') bins. Inner error bars show the statistical uncertainty and the total error bars represent statistical and systematic 
uncertainties added in quadrature. The parameterization of the curves is given by Eq. (|66|) and their parameters are given in 
Tab. m Central lines are calculated with the fitted values of parameters, while the dashed lines correspond to one standard 
deviation of the curve parameter. Except for the second —t' bin, the data points are shifted for better visibility. 




[GeV^] [GeV^] 

Fig. 5. The dependence of Re(Tii/Too) (left panel) and Im(Tii/Too) (right panel) for proton and deuteron data. Points 
show the amplitude ratios from Tabs. [2] and [3] after averaging over four —t' bins using Eqs. (|67|) and (|68|) . The parameterization 
of the curves obtained from combined proton and deuteron data is given by Eqs. (|66|) and (|69p while the parameters are given 
in Tab. [l] The meaning of the error bars and the explanation of the curves are the same as for Fig. |4] Here and hereafter 
proton data points are slightly displaced to the right for better visibility. As explained in the text, curves are shown as solid 
(dash-dotted) lines if the parameterization is based (not based) on pQCD predictions. 
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target 


ratio 


a, GeV 


Sa, GeV 




proton 


Re(rii/roo) 


1.173 


±0.048 


0.63 


deuteron 


Re(Tii/Too) 


1.106 


±0.035 


1.04 


proton+dcuteron 


Re(Tii/roo) 


1.109 


±0.026 


0.78 


target 


ratio 


b, GeV"' 


Sb, GeV"' 




proton 


Im(Tn/Too) 


0.340 


±0.025 


0.54 


deuteron 


Im(Tii/ro,j) 


0.359 


±0.020 


0.71 


proton+deuteron 


Im(Tii/roo) 


0.348 


±0.017 


0.75 



Table 4. The dependence of Re(Tii/Too) and Im(rii/roo) for proton, deuteron and combined data sets parameterized as 
given by Eqs. (|66|l and (|69|) . The values of parameters with their total uncertainties are presented. The last column shows the 
value of per degree of freedom. 



statistical uncertainty STstat where Om is the statistical 
uncertainty of Tm- 

In the left panel of Fig. [SI the average values of the 
ratio Re(tii) in each of the bins for both the proton 
and deuteron targets are shown. As the results for the two 
targets are found to be compatible, a fit to the combined 
proton and deuteron data is performed and the result ex- 
tracted for a in Eq. is also presented in Tab. S] Rea- 
sonable values of per degree of freedom, /Ndf ~ 1, 
are obtained. A comparison of the average values of ra- 
tio Re(tii) across —t' bins and the curve calculated with 
Eq. (|66p using the value of the parameter a obtained from 
the combined proton and deuteron data is shown in the 
left panel of Fig. O The dependence of Re(tii) is found 
to be in a good agreement with the asymptotic behavior 
expected from pQCD. 

As shown in the right panel of Fig. [51 the amplitude 
ratio Im(tii) rises with . This dependence of Im(tii) 
disagrees with the prediction [301(31] given by Eq. ([Fl]) . 
whereas the dependence of the real part of the same ratio 
^11 agrees with predictions based on the same formula. A 
fit of the function in Eq. (|66|) to the experimental data 
does not provide reasonable values of pcr degree of 
freedom: x^ /Ndf = 2.46 for the proton data and x^ /Ndf = 
4.25 for the deuteron data. Instead, a fit function of the 
type 

Im(tn) = bQ (69) 

gives reasonable results as can be seen from Tab. [H where 
fit results for the proton, deuteron and combined datascts 
are presented. The right panel of Fig. [SI shows the curve 
calculated with the parameter b obtained for the combined 
proton and deuteron data compared with the —t' bin av- 
eraged ratio Im(iii) for both proton and deuteron data. 
The disagreement with the prediction of Eq. (|61[) may be 
due to the fact that the range of the HERMES experi- 
ment is not in the asymptotic region. This is in agreement 
with the conjecture of Ref. [HI] that the hard scale for 
vector-meson electroproduction is Qv given by Eq. (|65p 
which is smaller than Q. Another possible explanation can 
be found in the discussion later in this paper of potential 
final-state interactions between the struck parton and the 
target remnant. 



The phase difference between the helicity ampli- 
tudes Til and Tqo is given by the relation 

tanJii = Im(tii)/Re(tii) . (70) 

It increases with as shown in the left panel of Fig. [H] 
Also shown is a fit to the dependence of 5ii calcu- 
lated with the functional form tan (5ii = bQ^ /a deduced 
from Eqs. (|551) and ([M]) . The result of the fit is Sn = 
(31.5 ± 1.4) degrees at (Q^) ^ 1.95 GeV^, obtained with 
the parameters a and b given in Tab. [4] for the com- 
bined proton and deuteron data. It is consistent within 
one standard deviation with the published result |18| ob- 
tained with the SDME method using the same proton and 
deuteron data as in the present analysis. A large phase 
difference Sn « 20 degrees was measured by the HI col- 
laboration for exclusive and (j) meson electroproduc- 
tion at 2.5 GcV2 < Q2 < 60 GcV^ [12]. A large value of 
(5ii = (33.0 ± 7.4) degrees was measured also in </> meson 
production at HERMES |ST] . As two-gluon exchange dom- 
inates in (j) meson production at HERMES kinematics, the 
measured large value of Sn for p'^ cannot be attributed 
solely to the quark-antiquark exchange essential for me- 
son production. This value of Sn is in clear disagreement 
with a calculation performed using a GPD-based ap- 
proach in pQCD that predicts a very small value for Sn. 
At present, there exists no model capable of explaining 
the value and dependence of 5ii. 

In existing pQCD calculations for exclusive vector- 
meson electroproduction, only two-gluon and/or quark- 
antiquark exchange have been taken into account. It can 
be argued |S3] that even in inclusive deep-inelastic lepton- 
nucleon scattering, the final-state interaction of the struck 
quark alters the cross section due to multi-gluon exchanges 
with the target remnant. This is at variance with the tra- 
ditional understanding obtained from inclusive DIS. It is 
shown |53] that a summation of multi-gluon-exchange am- 
plitudes results in an eikonal-like correction. For vector- 
meson production, if rescattering on the nucleon occured 
also for the quark-antiquark pair that transforms to the 
final vector meson, the eikonal-like correction might be 
responsible for the measured large phases observed in the 
present work. 

Figure [71 shows the t' dependence of the real and imag- 
inary parts of the ratio tn. Since Re(<ii) and Im(tii) de- 
pend on according to Eqs. ([M]) and ([5^ . they are 
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Fig. 6. The dependence of the phase difference &\\ (left panel) and 5oi (right panel, see Sec. l7.4p between the amplitudes T\\ 
and Toi, respectively, and Too obtained for proton and deuteron data. Points show the phase differences 5ii and (5oi calculated 
from ratios of amplitudes given in Tabs. [2] and [3] after averaging over —i! bins. Inner error bars show the statistical uncertainty 
and the outer ones show the statistical and systematic uncertainties added in quadrature. The fitted parameterization is given 
by Eqs. ()70|) and (|78l) respectively for i\\ and (5oi. The parameters of the curves are given in Tabs. |4] and [6] for combined proton 
and deuteron data. The central lines are calculated with the fitted values of the parameters, while the dashed lines correspond 
to one standard deviation in the uncertainty of the curve parameter. 
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Fig. 7. The t' dependence of Q ■ Re(Tii/Too) (left panel) and Im(rii/Too)/Q (right panel) for proton and deuteron data. 
Points show the amplitude ratios from Tabs. [5] and |3] after averaging over four bins using Eqs. (|67|l and (|68|l . The straight 
lines in the left and right panel show the value of a and fe, respectively, from Eqs. H66|) and (|69|l while the parameters a and h 
are given in Tab. |4] The meaning of the error bars and the explanation of the curves are the same as for Fig. |4l 
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Fig. 8. The dependences of |L*'ii/Too| on and t' for proton and deuteron data. The points show the amphtude ratios given 
in Tabs. [2] and [3] after averaging over —t' (Q^) bins in the left (right) panel. The inner error bars show the statistical uncertainty 
and the outer ones show the statistical and systematic uncertainties added in quadrature. The results fitting the combined data 
set with a constant (central line), |(7ii/Too| — g, are given in Tab.[Sl The dashed lines correspond to one standard deviation in 
the total uncertainty. 



target 


ratio 


9 


Sg 




proton 


1(711 /Tool 


0.400 


±0.020 


0.60 


deuteron 


It/ii/Tool 


0.383 


±0.017 


0.40 


proton+dcutcron 


\Uii/Tm\ 


0.390 


±0.013 


0.49 



Table 5. Results of fitting the ratio |i7ii/Too| to a constant for proton, deuteron and combined data sets. The values of 
parameters with their total uncertainties are presented. The last column shows the value of per degree of freedom. 



shown in Fig. [7| multiplied or divided by Q, respectively. 
No noticeable t' dependence is observed for Re(tii) and 
Im(tii). Since the differential cross section of the pro- 
cess in Eq. ^ for high energies and small \t'\ is usu- 
ally described by an exponential factor exp{/3t'}, the he- 
licity amplitudes should have exponential factors Tqq cx 
exp{l3Lt'/2} and Tn cx exp{i3Tt'/2}. The absence of a t' 
dependence of the ratio Tn/Too means that the slope pa- 
rameters /3l and /3t for the amplitudes Too Stiid Tn are 
close to each other. For very small \t'\, it is reasonable to 
use the linear approximation 



Rc(rii/Too)= -cxp{--zi/3it'} 

^(l±izi/3i|t'|), (71) 

Im(Tn/roo)- 6Q cxp{-iz\/32t'} 

bQ{l + ^Al32\t'\). (72) 



The proton results are A^i = (-1.02 ± 0.85) GeV"^ ^nd 
Z\/32 = (-0.91 ± 2.00) GeV-2, while the fit for the deute- 
ron data gives = (0.58 ± 0.80) GeV"^ and Z\^2 = 
(-1.96 ± 1.58) GeV-2. The results of the fits show that 
all four numbers are consistent with one another. We now 
assume that, within experimental accuracy, the slope pa- 
rameters for the real and imaginary parts of the ratio 
coincide across both target types. In this case we have 
APi « A/32 ~ Pl — Pt- Combining these four numbers 
making use of Eqs. ([ST]) and we get an estimate for 
Pl- = (-0.4 ± 0.5) GeV-2. This result on Pl - Pt 
is in agreement with the prediction published in Ref . |54) , 
which ranges from -0.7 GeV'^ at = 0.8 GeV^ to -0.4 
GcY-^ at g2 5 GcV^. 



7.3 Kinematic Dependence of |C/ii/Too| 

The unnatural-parity-exchange amplitude, Uu, describes 
the transition from a transversely polarized photon to a 
transversely polarized p meson 

(7t ^- Pt)- At large W 
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0.97 


deuteron 


Im(Toi/Too) 


0.161 


±0.090 








1.20 


proton+deuteron 


Ini(Toi/Too) 


0.221 


±0.069 








1.10 


target 


ratio 


/i GeV"' 


Sfi GeV-i 


/2 GeV"^ 


Sf2 GeV-3 


Pc 




proton+deuteron 


Im(Toi/Too) 


0.653 


±0.132 


-0.285 


±0.065 


-0.903 


0.66 



Table 6. The kinematic dependence of Re(roi/roo) and Im(Toi/Too) for the proton, deuteron and combined data sets 
parameterized as given by Eqs. (|74I - I76|I . The values of parameters with their total uncertainties are presented. The last column 
shows the value of P^r degree of freedom and pc is the correlation parameter. 



and , this transition should be suppressed by a fac- 
tor of My/Q compared to the dominant amplitude Tpo 
as explained in Sec. 14.2.21 The UPE contributions to the 
amplitude may be sizable at intermediate energies 
typical for HERMES. 

The ratio |uii| versus and t' is presented in Fig. |51 
The value of |J7ii| is found to be smaller than |Too| by 
a factor of approximately 2.5. No kinematic dependences 
of the ratio |uii| are seen and therefore it is fitted to a 
constant 

\uii\=g. (73) 

The results of the fit to proton and deuteron combined 
data are given in Tab. [51 and the constant g is shown 
in Fig. [5] by straight lines. Fits with the same function 
a/Q as for Re(tii) (see Eq. corresponding to the 

behavior expected in pQCD, give values for pcr degree 
of freedom 2.03, 2.67, and 2.25 for the proton, deuteron, 
and combined data, respectively. This disagreement may 
reflect the fact that the HERMES region is far from 
the asymptotic one. 

From a study of soft hadron scattering, it is known j35l 
[551155] that the most important contribution to UPE am- 
plitudes at intermediate energies comes from pion excha- 
nge. The amplitude of one-pion-exchange in the t-channel 
contains the propagator l/(t — mj), which become large 
for small values of — t that approach the pole at the un- 
physical value of — t = — ~ —0.018 GeV^. The mean 
value of —t' for the first —t' bin is 0.019 GeV^, and hence 
might appear to be small enough to approach the vicinity 
of the pole. However, as shown in the right panel of Fig. [8l 
no evidence for a pion-pole-like dependence can be seen 
in the data within the statistical precision of the measure- 
ment. Such a dependence could be weakened by various 
effects: 

i) even in the first bin in —t', the value of {—t) equal to 
0.06 GeV^ is substantial, and hence not close enough to 
the pion pole at unphysical positive t, 

ii) the amplitude Tqq has a strong exponential dependence 
cx exp{l3Lt'/2} with /3l « 7 GcV^^ and decreases 
with |t'|, hence the ratio Uh/Tqq is flatter than Uu itself. 



iii) in addition to one-pion exchange, other exchange pro- 
cesses can contribute to the amplitude Uu. 

Using the amplitude method, the signal of unnatural- 
parity exchange has a significance of more than 20 stan- 
dard deviations in the total experimental uncertainty sep- 
arately for each of the proton and deuteron data sets (see 
Tab. [51 and Fig. [5]). In contrast, the existence of UPE was 
established [TS] with a significance of 3 standard devia- 
tions for the combined proton and deuteron data in the 
analysis using the SDME method. 

7.4 Kinematic dependence of Toi/Too 

The amplitude Tqi = ^oiii describing the transition 7^ — >■ 
is expected to be the largest SCHC-violating ampli- 
tude. In accordance with the asymptotic formula (|62p . the 
parameterization 

Im(ioi) = (74) 

is used. A fit to both proton and deuteron data using this 
parameterization gives acceptable values, as seen in 
Tab. [51 Proton and deuteron results are compatible within 
one standard deviation in the total uncertainty, although 
neither measurement is particularly precise as shown in 
Tab. [51 and the right panel of Fig. [SI in which the data are 
multiplied by Q in order to demonstrate the \J —t' depen- 
dence. The fit of the combined proton and deuteron data 
set yields a value of the parameter / that is three standard 
deviations from zero with /Ndf ~ 0.66 (see Tab.[ni). We 
note that these results do not necessarily confirm the va- 
lidity of Eq. ((52) although the data do not contradict the 
pQCD prediction. 

The dependence of the amplitude ratio Im(ioi ) for 
the proton and deuteron data is shown in the right panel of 
Fig-Hni As shown in the figure, the point for the deuteron 
in the largest bin is slightly negative which favours a 
fit using the equation 

Im(<oi) = \/^(/i+/2g') (75) 
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Fig. 9. The t' dependence of Re(roi/roo) (left panel) and Q ■ Im(roi/roo) (right panel) for proton and deuteron data. Points 
show amplitude ratios from Tabs. [2] and [3] after averaging over bins. Inner error bars show the statistical uncertainty and 
the outer bars show statistical and systematic uncertainties added in quadrature. The parameterization is given by Eqs. (|74|) 
and (|76|l . The parameters of the curves are given in Tab. [S] for combined proton and deuteron data. Central lines are calculated 
with fitted values of parameters, while the dashed lines correspond to one standard deviation of the curve parameter. 
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Fig. 10. The dependence of Re (To i /Too ) / V ~ i ' (left panel) and Im(roi /Too)/ ^/—t' (right panel) for proton and deuteron 
data. Points show amplitude ratios from Tabs. [2] and [3] after averaging over —t' bins. Inner error bars show the statistical 
uncertainty and the outer bars show statistical and systematic uncertainties added in quadrature. The parameterization is given 
by Eqs. (|75p and (|76|l . The parameters of the curves are given in Tab.[S]for combined proton and deuteron data. Central lines are 
calculated with fitted values of parameters, while the dashed lines correspond to one standard deviation of the curve parameter. 
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over that using parameterization ([71]). The resuh of the fit 
with function ((75)) is shown in the right panel of Fig. [TOl 
The parameters /i and /2 are strongly correlated. The 
correlation parameter pc is presented in Tab. |5] and is 
taken into account in the calculation of the uncertainty 
of Im(ioi) using Eq. ([75|. 

The dependence of the ratio Re(<oi)/'\/— is pre- 
sented in the left panel of Fig. [TO] The quantity Re(toi) 
l\J —t' does not decrease with and can be well de- 
scribed by a constant which is also shown in the figure. 
Using for Re(ioi) the functional form of Eq. (fM)) yields 
/Ndf values of 1.08, 1.55, and 1.31 for separate fits to 
the proton, deuteron and combined data sets. Using in- 
stead the simpler parameterization 

Re(ioi) = c^f^ (76) 

decreases the value of by a factor of approximately 
two, indicating a better description of the data. The re- 
sults of this fit are shown in Tab. [5] and the left panels of 
Figs, ini and [TUl The small values of /Ndf as shown in 
Tab. [5] indicate that our systematic uncertainty might be 
overestimated. 

The phase difference Jqi between the amplitudes Tgi 
and Too can be calculated using the formula 

tan,5oi =Im(iio)/Re(tio). (77) 

The results for the dependence of Jqi after averaging 
over —t' bins are shown in the right panel of Fig. |6l The 
phase difference (5oi is non-zero and decreases with Q^. It 
is slightly negative for the deuteron data at the largest 
value of Q^. Using instead the results from the fits to 
formulas ([75l [75]) given in Tab. [HI the phase difference is 
obtained as 

tanJoi = (/i + /2Q')/c, (78) 

which is shown as a curve in the right panel of Fig. [HI As 
can be seen, the fit function ([75]) is able to describes the 
dependence of the phase difference (5oi- 



7.5 Small Amplitude Ratios Tio/Too and Ti i/Too 

According to the hierarchy given in Eq. ([5T]) . Tio and Ti-i 
are the smallest amplitudes that give linear contributions 
to the numerators of SDMEs in the case of an unpolar- 
ized target. The ratio tio is expected to be proportional 
to ^/—t' in accordance with Eq. ([60]) . while is ex- 
pected to be proportional to t' . These ratios are supposed 
to depend on in the asymptotic region according to 
Eqs. ([63]) and ([64])- Figure [TT] shows Re(tio)/\/^ and 
Im(tio)/\/--F versus for both proton and deuteron 
data. As shown in Fig. [Til the values of Re(iio) / \/— i' and 
Im(tio)/ yf—t' on the proton are compatible with zero. The 
results of the fit of Re(iio) with the function 

Re(iio) = (79) 



and Im(iio) with the same function are presented in Tab. [7] 
The ratio Re(<io )/%/—<' for the deuteron is slightly pos- 
itive (except the first point at the smallest value of Q^) 
and increases with . The fit of the deuteron data with 
the function 

Re(tio) = sQ^^T^' (80) 

provides the positive result which is presented in Tab. [7] 
As shown in Fig.[TTl the quantity Im(tio)/ \f—V is negative 
for the three points at the largest values of . The fit of 
Im(tio) / \J using the same function (|80p gives the neg- 
ative result presented in Tab. [7] This behavior contradicts 
the pQCD prediction given by Eq. (|63p . and indicates that 
in the case of the amplitude T\q the HERMES kinematic 
range may be far from the large-Q asymptotic region. 

We recall that the amplitudes Too and Tio vanish in 
the small Q-limit (Q^ — )■ 0) because the real photon has 
no longitudinal polarization. If the behavior at fixed i! and 
small were described by 

Too « Q , (81) 
Tio « Q^ (82) 

the behavior of the ratio t\Q = Tiq/Tqq would be just as 
assumed by Eq. ([50]) and presented in Fig. [TT] 

The ratios Re(ti_i)/(— t') and Im(ii_i)/(— t') are pre- 
sented versus in Fig. [T^ Both ratios appear to be com- 
patible with zero within 2.5 a of the total uncertainty both 
for the proton and deuteron target and only the combined 
data provides a non-zero signal for Im(ii_i). The fit of 
Re(ti_i) with the function 

Re(ti_i) = i-t')^ (83) 

corresponds to the second term in Eq. ([M]) . The result, 
presented in Tab. [S] for the combined proton and deuteron 
data, shows a signal of 2(t significance with respect to the 
total uncertainty. The result of the fit of Im(ii_i) with 
the same function ([55]) is also presented in Tab. [5] The pa- 
rameter h deviates from zero by more than three standard 
deviations of the total uncertainty. A fit of the combined 
proton and deuteron data with both terms in Eq. (|64[) 
provides values of Ci and C2 that are consistent with zero 
within 1.5 cr of the total uncertainty for both Re(ti_i) 
and Im(ii_i). This means that the statistical precision of 
the HERMES data does not allow a reliable verification of 
Eq. ([M]) or an extraction of the gluon transversity GPD. 

7.6 Comparison of HERMES and HI Results 

Here we compare the results of the present work with the 
analysis of the p^-meson production data on the proton 
by the HI collaboration [5D], obtained for the CM en- 
ergy range 36 GeV < W < 180 GeV, photon virtuality 
2.5 GeV2 < Q2 < 60 GeV^, and -t' < 3 GeV^. For the 
comparison of the t' dependence, the results at the mean 
value of equal to 3.3 GeV^ are chosen which arc the 
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target 


ratio 


r, GeV-i 


Sr, GeV^^ 




proton 


Ro(rio/Too) 


-0.012 


±0.030 


1.07 


proton 


im(rio/roo) 


0.019 


±0.061 


0.58 


target 


ratio 


s, GeV"^ 


5s, GcV-^ 




deuteron 


Re(rio/Too) 


0.045 


±0.015 


0.32 


deuteron 


im(rio/roo) 


-0.109 


±0.021 


1.02 



Table 7. The kinematic dependences of Re(rio/roo) and Im(rio/Too) for proton and deuteron parameterized as given by 
Eqs. (|79p and (|80|l . The values of parameters with their total uncertainties are presented. The last column shows the value of 
per degree of freedom. 



target 


ratio 


h, GeV"' 


5h, GeV~^ 




proton+deuteron 


Re(Ti_i/roo) 


-0.156 


±0.059 


0.79 


proton±deuteron 


im(ri_i/roi,) 


-0.418 


±0.126 


0.54 



Table 8. The kinematic dependence of Re(ri_i/Too) and Im(ri_i/Too) for combined proton and deuteron data parameterized 
as given by Eq. (|83|l . The values of parameters with their total uncertainties are presented. The last column shows the value of 
per degree of freedom. 




Fig. 11. The dependence of Re(Tio/roo)/\/~i' (left panel) and Im(rio /Too)/ \/—t' (right panel) for proton and deuteron 
data. Points show amplitude ratios from Tabs. [5] and O after averaging over —t' bins. Inner error bars show the statistical 
uncertainty and the outer ones indicate statistical and systematic uncertainties added in quadrature. The parameterization of 
the curves describing the deuteron data both on Re(rio/roo) and Im(rio/Too) is given by Eq. (|80[l and their parameters are 
given in Tab. [T] Central lines are calculated with the fitted values of parameters, while the dashed lines correspond to one 
standard deviation of the curve parameter. 



closest ones to the region of the HERMES data. In 
the HI data analysis presented in Ref. [20], the imaginary 
parts of the amplitude ratios were not extracted. The addi- 
tional hypothesis was used in Ref. [20] that for any ratio 
of the amplitudes, the approximate relation |tAvA^P = 
[Re(fAvA.^)]^ is valid within the experimental accuracy. 
The comparison of the HERMES and HI results is pre- 
sented in Fig. [131 As seen from the figure the HERMES 



and HI results agree within their total uncertainties. No 
strong dependence on W is observed for the amplitude 
ratios Rc(tii) and Re(ioi). Considering the differences 
in Q2 between HERMES ((Q^) = 2.0 GeV^) and HI 
((Q^) = 3.3 GeV^), no dependence can be seen for 
Rc(toi) at small even for < 3.3 GeV^ (see also the 
discussion in Sec. 17. 4p . The ratios Re(tio) and Rc(ii_i) 
obtained by HI [5D] are compatible with zero at = 3.1 
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Fig. 12. The dependence of Re(Ti_i/Too)/(-t') (left panel) and Im(Ti_i/Too)/(-i') (right panel) for proton and deuteron 
data. Points show amplitude ratios from Tabs. [2] and [3] after averaging over —t' bins. Inner error bars show the statistical 
uncertainty and the outer ones indicate statistical and systematic uncertainties added in quadrature. The parameterization of 
the curves describing the combined proton and deuteron data both on Re(Ti_i/Too) and Im(ri_i/roo) is given by Eq. ()83p and 
their parameters are given in Tab. [S] Central lines are calculated with the fitted values of parameters, while the dashed lines 
correspond to one standard deviation of the curve parameter. 
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Fig. 13. The kinematic dependences of Re(rii/roo) and Re(roi/roo) for proton and deuteron data. Points of HERMES are 
the same as in Figs. [5] and [9l The HI results are from Ref. [20]. Inner error bars show the statistical uncertainty and the outer 
ones indicate statistical and systematic uncertainties added in quadrature. 



GeV^ within the experimental accuracy. This is consistent with the HERMES results for the proton shown in Figs.fTTI 

andm 
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8 Conclusions 

Exclusive p'^-meson electroproduction is studied in the 
HERMES experiment, using a 27.6 GeV longitudinally po- 
larized electron/positron beam and unpolarized hydrogen 
and deuterium targets in the kinematic region 0.5 GeV^ < 
g2 < 7.0 GeV2, 3.0 GeV < < 6.3 GeV, and -t' < 
0.4 GeV^. These data are used to determine the real and 
imaginary parts of the ratios Tn/Too, Toi/Too, Tiq/Tqo, 
Ti_i/Too, and |;7ii/Too| for 16 bins in and -t'. System- 
atic uncertainties due to the background contribution, un- 
certainties in the Monte Carlo input parameters, and the 
uncertainty of the applied amplitude method are given. 
Except for Tiq/Tqq, the amplitude ratios for the proton are 
compatible with those for the deuteron. The extracted am- 
plitude ratios Tu/Tqo, \Uii/Tqq\ and Tqi/Tqq are found to 
be sizable. The ratios Tiq/Tqq and Ti-i/Tqo for the proton 
are found to be compatible with zero within experimental 
uncertainties. The ratio Ti_i/Too for the deuteron is also 
zero within experimental accuracy, while Re(Tio/T'oo) is 
slightly positive and Im(Tio/T'oo) is slightly negative, ex- 
cept for the bin with the smallest values. 

The SDMEs calculated in terms of these helicity am- 
plitude ratios generally agree with the results published 
in Ref. jl8j . The amplitude method was shown to provide 
more accurate polarized SDMEs than the previous analy- 
sis using the SDME method. 

The statistical precision available in this analysis per- 
mits the parameterization of the kinematic dependences 
of amplitude ratios and the extraction of the phase dif- 
ference between various helicity amplitudes. The real part 
of Tu/Tqq is found to follow the asymptotic 1/Q behav- 
ior predicted by pQCD [5IJll31j . The imaginary part of 
Tii/Tqq grows with Q^, in contradiction to the large-Q^ 
asymptotic behavior expected from pQCD. The phase dif- 
ference Sii between the amplitudes Tn and Too grows 
with and has a mean value of about 30 degrees in the 
HERMES kinematic region. This is in agreement with the 
published result of the SDME method [TH] and in contra- 
diction to calculations |lllll2l[T3l l30 l [3T ] based on pQCD. 
For the first time, the dependence of Su observed in 
Ref. [IH] is shown to be related to the increase with of 
the imaginary part of the ratio of the helicity amplitudes 
Tii/Tqq. 

The behavior of Im(Toi/Too) is found to be in agree- 
ment with the asymptotic pQCD behavior V—t'/Q, while 
the extracted value of Re(Toi/Too) is likely to be in dis- 
agreement with the pQCD prediction. The data indicate 
non-zero values of the phase difference Sqi between the 
amplitudes Tqi and Too and the decrease of Sqi with Q^. 

The ratio |?7ii/Too| is found to be constant in the HER- 
MES kinematic region, in disagreement with the asymp- 
totic pQCD behavior at large Q^. No pion-pole-like behav- 
ior is observed at small \t'\. The UPE signal is seen with 
a very high significance for both proton and deuteron data, 
confirming the existence of unnatural-parity exchange con- 
tributions with a higher precision than that obtained with 
the SDME method [18]. 

The dependence of the amplitude ratio Re(Tii/Too) 
and t' dependence of the amplitude ratio Re(Toi/Too) are 



also compared to those extracted by the HI collabora- 
tion at the center-of-mass energy range 36 GeV < W < 
180 GeV, photon virtuality 2.5 GeV^ < < 60 GeV^, 
and —t' < 3 GeV^. No strong dependence of the ampli- 
tude ratios Re(Tii/Too) and Re(Toi/Too) on W, and no 
dependence of Re(ioi) at small \t'\ are observed. The 
HERMES and HI results agree within their total uncer- 
tainties. 
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Appendix 

A Systematic Uncertainty of the Amplitude 
Method for an Unpolarized Target 

The systematic uncertainty arising from the neglect of am- 
plitudes that make small contributions to the angular dis- 
tribution applies if only data on unpolarized targets are 
considered. If the set of observables were enlarged and 
data on targets with transverse and longitudinal polariza- 
tions were added, then all the amplitude ratios could be 
extracted and there would be no contributory systematic 
uncertainty of the amplitude method. 

The SDME analysis ^ of the HERMES data has 
shown the absence of any signal of the UPE amplitudes 
Uqi, C/io, and violating the SCHC approximation. 

The contribution of the greatest UPE amplitude Un to 
any SDME changes it by a value which is less or about 
one standard deviation of the statistical uncertainty. We 
neglect the contributions of Uqi, Uiq, and C/i_i to the 
SDMEs and do not consider that as a possible source of 
systematic uncertainty. The terms with the amplitudes 
f^oij Uio, and f/i_i in all the formulas considered in this 
appendix are also ignored. As explained in Sec. 14.2.11 the 
fractional contribution to SDMEs of NPE helicity-flip am- 
plitudes are suppressed by . The systematic uncertainty 
of the extracted amplitude ratios due to the neglect of the 
NPE nucleon helicity-flip amplitudes is the only uncer- 
tainty of the amplitude method which is considered below. 

The true ratio of the NPE amplitude without helicity 
flip T^yi^^i to TqIqi is denoted by tx^x^. The true ratio 
of the nucleon-helicity-flip amplitude T^^^ = T^^^^_i-)^^i 
to TqIqi is denoted by tl^x.,^ in particular, tgo = T^-^o^ 

/T,.,., while |«nP ^ (It^iiiP + |C/i-iiiP)/|T},ioiP- 
In the following, we estimate the effect of neglecting small 
contributions of the NPE nucleon spin-flip amplitudes, 
-^VA ' which bias the fitted amplitude ratios by Stx^-x^ 
and respectively. In order to avoid misunderstand- 

ings, we note that the notation tx^x^ and |uii| in the 
main text was used for the fitted amplitude ratios which 
corresponds now to txyX-, + ^tx^X-, and |mii| -I- i5|ziii|, re- 
spectively. 

As shown in Eqs. (|24II26|) . any SDME is multiplied by 
a function of the angles 9, <P, (p and these functions are lin- 
early independent. This means that every SDME (rather 



26 



than combinations of them) is determined in the fit indi- 
vidually. Now we assume that the true amplitude ratios 
are known, so that the true values of all SDMEs can be 
calculated. 

The exact expression for Tqq, which was presented in 
Ref. [IH] and rewritten here in terms of the ratios of the 
helicity amplitudes tx^-\^, t\^,^ and is 



eii+mo\') + \toi\' + \tU' 



'00 



where 



e 1 + It; 



' |2 + 2|tinl^ + 2 



(Al) 



\tlo\') 



+ |iiir + |tnr + |ioi|^ + |tgil 

2 



+ \tl-l\' + \tU\'' 



"11 



(A2) 



It is more convenient technically to consider a logarithm of 
Tqq. The contribution to In(rQQ) of the small NPE nucleon- 
helicity-flip amplitudes Itg^p, jtfop, |tf jtg^p and |tf_iP 
linear in these small quantities, is 



s 12 



- 'ei\tt,o\' + 2\tl,\') + \tl,\' + \t^,,\^ + \tU\' 
I e(l + 2|iion + |<n|' + |toiP 



+ 1^1-1 
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(A3) 



The formulas for Tqq used in the fit of the angular 



distribution are 



'00 



e+|<oig + '^|toi|' 



(A4) 



where 

Ntot 



i + 2(^itior + <5|iioi 

+ (|tiiP + '5|inp) + (|<oi|' + '^|ioir 
+ (|ti-iP + b\ti-xf) + (l^^iiP + &Wxf) {kh) 

with (5|tAvAj^ = ^Av-A.^'^^AvA-, + ^Av^A., <^tAv A^ ^nd (5|uii|2 

= 2|uii|(5|uii|. 

The contribution to In(rQQ) linear in the small quanti- 
ties ^\t\^\ P and 5|uAv-A-yP is given by the relation 



<51n(?«4) 

2e(5|tio 
+ <5|ii-i 



'00 
' + <5|tllP 



e+|ioiP 
' ' %iiP 



-^ItoiP 



/ [e(l + 2|tiop) +|tiiP + |toi|' 



+ |ti-ip + |unp 



(A6) 



A comparison of Eqs. (|A1I) and (|A2|) with Eqs. (|A4|) and 



Av A^ 



(jA5|) shows that r{}g = when 8t\^\_^ = | = t 
0. Since the fit has to reproduce the angular distribution, 
the contribution of the nucleon-helicity-flip amplitudes to 
Tqq is compensated by the deviation of the obtained ampli- 
tude ratios from the true amplitude ratios. This assumes 
the validity of the relation 



A\T,{rl\)^6\^{l%\) (A7) 
which provides the linear relation between \t\^y^ P and 
b\txy\^ |2, 5\uxv\^ p. Considering all other SDMEs r^^^,^, 
use of the same method provides the set of relations for 

7/ = 1 - 3, 5-8 

^ln(rLA'.,)-'^ln(^,,A'„) (A8) 



which can be used to determine St\^-x^ and (5|wii|. 

In order to get an approximate solution of Eq. (|A7[) . 
the contributions of the s-channcl helicity violating am- 
plitudes Toi, Tio, and Ti_i {T^^, Tf^, and Tf.i) are ne- 
glected since they are small compared to Too s^nd Tn (Tqo 
and Tfi). Putting Eqs. ^M^l and ^M\i into Eq. (|SZ| yields 
the relation 



s |2 
00 1 



s 12 
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<^|tl 



S\uii 



e+ |<ii|2 + |uii| 



e+|iiiP + |uii| 



(A9) 



Considering in the same way the matrix element rl_i (or 



Im(r^_j^)) and using Eq. (|A8I) we get the relation 



lifiP 



elt; 
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l"lll^ 



6\t 



111 



l"ii|^ 
<5|«iiP 



111 



|Mii| 



e + t 
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The solution of the system of two equations 
(rnoll is 



S\t 
S\ui 



111 



1^11 
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(All) 
(A12) 
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Dividing Eq. (|IlT|) by \tii\^ and Eq. Pl2|) by 
solutions can be rewritten for the fractional systematic 
uncertainties of l^nP and |miiP in the form 



S\tn\' 
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S\uii\'^ ^S\uii\ 
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(A14) 



In order to obtain estimates for the real and imaginary 
parts of 6tii, we consider the SDMEs Re(rfo) and Re(rfo), 
respectively. The exact expression for Re(rfo) taken from 
Ref. [T5] and rewritten in terms of ^Av-a.^, ^va ' ^'^'^ l"iil 



Reirlo) - -^Re 



tn -t 



1-1 



+ ifii-tuMor 



2<io(<oi)* + 2i^o(^gi)* /Ntot (A15) 
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where Ntot is defined by Eq. Using Eq. ^ an Since \Rc{5tx^x^)\ < \5tx^,x^\ and |Im(5tAv-A^)l 

approximate equation analogous to Eqs. (jA9p and (jAlOp < |(5iAvA^|, the right hand side of relation (jA26p can be 

used both for the real and imaginary parts of Stx^x^ yield- 
ing Eqs. ([57)) and Relation ([5^ follows immediately 
from Eqs. (|M4|) and ((X24)) . 



is obtained: 



Rc[^iiftgo)* 
Re(tn) 

Re{6tn) 
Re(tii) " 



e+ KiiM + l"ii| 
<5|iiiP + J|?.nP 
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(A16) 



Substituting into this equation the solutions for p and 
given in Eqs. (jAlip and (|A12|) leads to: 

Re(<5iu) = RcKi(tgo)*] - \ej'Rc{tn). (A17) 

Considering Re(rfQ) analogously we get 

Im(Jtn) = ImKi(ig„)*] - |igopIm(in). (A18) 

Combining Re((5tii) and Im((5tii) into the complex num- 
ber 5^11 = Re((5tii) + i ■ Im((5iii), the equation 
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(A19) 



is obtained which is equivalent to Eqs. (jA17|) and (jA18|) . 
From Eq. (Pl9)) . it follows that 
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An analogous consideration of Tqq and 7'qq provides the 
expression 



^_toi^rj^_ ioo\riooY 

toi ^Tqi Tf)QJ \TqqJ 
while considering r\i and Im(ri^_]^) leads to: 
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Consideration oilm{rf_i) and Im(?-[_]^) leads to the rela- 
tion 



tio 

As shown in Ref. 
obey the inequality 
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the ratios of T 



A v Ao 



and Ta^,a^ 



\n,xjTx,xJ<VT/i2M), 
which gives the inequality 
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Inserting Eqs. (|A24|) and (|A25P into expressions (|A20I 
IA23P yields the formulas of interest for the systematic 
uncertainty of tAvA., 
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